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Abstract 

The dependence of Tc on hydrostatic (He-gas) pressure for superconducting 
MgB2 has been determined to 0.7 GPa. We find that Tc decreases linearly 
and reversibly under pressure at the rate dTc/dP ~ —1.11 it 0.02 K/GPa. 
These studies were carried out on the same sample used in earlier structural 
studies under He-gas pressure which yielded the bulk modulus B = 147.2 it 0.7 
GPa. The value of the logarithmic volume derivative of Tc is thus accurately 
determined, d\nTc/dlnV = +4.16 it 0.08, allowing quantitative comparison 
with theory. The present results support the emerging picture that MgB2 is a 
BCS superconductor with electron-phonon pairing interaction. 



The recent discovery of superconductivity in MgB2 at ~ 39 K has sparked 
worldwide a torrent of experimental and theoretical activity, reminiscent of the frenzy 
following the observation of superconductivity in La-Ba-Cu-0 at a comparable 
temperature almost 15 years ago. Replacing ^''B with ^^B results in a sizeable 
isotope shift to lower temperatures which points to BCS superconductivity. Other 
experiments, such as heat capacity[^, ^, photoemission spectroscopy 0, and inelastic 
neutron scattering!]^, ^ also support the picture that MgB2 is a phonon-mediated 
superconductor in the weak-to-moderate coupling regime. 

High pressure studies traditionally play an important role in superconductivity. 
A large magnitude of the pressure derivative dTd dP is a good indication that higher 
values of T^. may be obtained through chemical means. It is not widely appreciated, 
however, that the pressure dependence Tc(P), like the isotope effect, contains valuable 
information on the superconducting mechanism itself. For example, in simple-metal 
BCS superconductors, like Al, In, Sn, and Pb, Tc invariably decreases under pressure 
due to the reduced electron-phonon coupling from lattice stiffening 0. More gener- 
ally, an accurate determination of the dependence of both and the lattice param- 
eters on pressure yields the functional dependence = Tc[a{P),b{P), c{P)] which 



provides a critical test of theoretical models. Hirsch[0 and Hirsch and MarsigliofTT 



have applied a theory of hole superconductivity to MgB2 and predicted that for an 
optimally doped sample Tc should increase with pressure, in contrast to the expected 
decrease in Tc from lattice stiffening. 

Precise structural data on MgB2 at room temperature (RT) have recently been 
obtained by Jorgensen et a/.|]T^ for hydrostatic pressures to 0.6 GPa in a He-gas 
neutron diffraction facility which yield the anisotropic compressibilities dlna/dP = 
-1.87 X 10-3 GFa-\ dhih/dP = -3.07 x lO'^ GFa~\ and the bulk modulus B = 
147.2 ±0.7 GPa; the compressibility along the c axis is thus significantly (64%) larger 
than that along the a axis. The binding within the boron layers is evidently much 
stronger than between the layers. These results are in reasonable agreement with 
electronic structure calculations by Loa and Syassen[0 . 



Recent synchrotron x-ray diffraction studies at RT in a diamond-anvil-cell (DAC) 



to much higher pressures (8 - 12 GPa) using dense He[|14[ or methanol-ethanol[T 
as hydrostatic pressure media yield the bulk moduli B = 155 ± 10 GPa and 151 ± 5 
GPa, respectively, in agreement, within experimental error, with the He-gas study [p!2 



A further DAC study|jT6|] to 7 GPa with silicon oil as pressure medium gives the 
significantly smaller value B = 120 ±5 GPa, even though the reported 2.9% decrease 
in unit cell volume upon applying 6.5 GPa pressure is less than the 4% decrease 
found in the other two DAC studies [p!4|, |15[. The strong deviation of the "silicon 



oil" data|]l6[ may arise from difficulties in extrapolating the data to zero pressure to 
obtain the bulk modulus and/or from shear stresses arising from the solidification of 
the silicon oil under pressure. The relatively large compressibility anisotropy in MgB2 
mandates the use of hydrostatic pressure in quantitative studies, since shear stresses 
applied by non-hydrostatic pressure media to an elastically anisotropic sample can 
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lead to erroneous results. 

Several studies of the dependence of Tc on pressure have appeared for MgB2. In 
the first experiment utilizing a fluid pressure medium (Fluorinert), Lorenz et al.\lT\ 
report that Tc decreases linearly with pressure to 1.8 GPa at the rate dTc/dP ~ —1.6 



K/GPa. On the other hand, using the same pressure medium, Saito et al. p8|] 
find the more rapid decrease dT^/dP ~ —2.0 K/GPa. Both groups cite their results 
to argue that MgB2 is a BCS phonon- mediated superconductor [^, |TB|, as is also 
argued by Loa and Syassen|l^. In further experiments to 25 GPa utilizing the solid 
pressure medium steatite, Monteverde et al. |]TU| find that Tc decreases under pressure 
at differing initial rates (-0.35 K/GPa to -0.8 K/GPa), Tc{P) showing a quadratic 
behavior for two of the four samples studied; each of the four samples was prepared 
using a different synthesis procedure. 

In the present experiment, Tc{P) is determined to 0.7 GPa in an ac susceptibility 
measurement using a He-gas pressure system; the helium pressure medium remains 
fluid at Tc ~ 39 K up to 0.5 GPa and thus applies true hydrostatic pressure (no 
shear stresses) to the sample. We find that Tc decreases linearly and reversibly with 
pressure at the rate dTc/dP ~ —1.11 ± 0.02 K/GPa. Implications for the nature of 
the superconducting state are discussed. 

The powder sample for this study was taken from the same mother sample used 
in parallel neutron diffraction studies [|12|. It is made using isotopically-enriched ^^B 
(Eagle Picher, 98.46 atomic % enrichment). A mixture of ^^B powder (less than 200 
mesh particle size) and chunks of Mg metal were reacted for 1.5 hours in a capped 
BN crucible at 800°C under an argon atmosphere of 50 bar. As seen in Fig. 1, the 
resulting sample displays sharp superconducting transitions in the ac susceptibility 
with full shielding and an onset temperature at ambient pressure Tc(0) ~ 39.25 K. 
Both x-ray and neutron diffraction data show the sample to be single phase with the 
AlB2-type structure. 

The present high pressure studies were carried out using a He-gas high-pressure 
system (Harwood) to 1.4 GPa; the pressure is determined by a calibrated manganin 
gauge in the compressor system at ambient temperature. The superconducting 
transition of the 8.12 mg MgB2 powder sample is measured by the ac susceptibility 
technique using a miniature primary/secondary coil system located inside the 7 mm 
I.D. bore of the pressure cell. A small Pb sphere with 1.76 mm dia. (38.58 mg) is also 
inserted in the coil system for susceptibility calibration purposes; for selected data 
the superconducting transition temperature of this Pb sphere serves as an internal 
manometer |]2D| to check the pressure indicated by the manganin gauge. The CuBe 
pressure cell (Unipress), which is connected to the compressor system by a 3 mm 
O.D. X 0.3 mm I.D. CuBe capillary tube, is inserted into a two-stage closed-cycle 
refrigerator (Leybold) operating in the temperature range 2 - 320 K. The pressure 
can be changed at any temperature above the melting curve Tm of the helium pressure 
medium (for example, T„ ^ 13.6 K at 0.1 GPa and ~ 38.6 K at 0.50 GPa[|Tl). For 



pressures above 0.5 GPa, lies above T^; the slight pressure drop (few 0.01 GPa's) 
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on cooling from to is estimated using the isochores of HeEl]]. All pressures 



are determined at the temperature Tc ~ 39 K. Further details of the experimental 
techniques are given elsewhere p^] . 

In Fig. 1 are shown representative examples of the superconducting transition in 
the ac susceptibility at both ambient and high pressure. With increasing pressure 
the narrow transition is seen to shift bodily to lower temperatures, allowing a deter- 
mination of the pressure-induced shift in Tc to within ± 10 mK. Remarkably, close 
inspection of the data for 0.50 GPa reveals a slight shift in the transition curve near 
its midpoint, accurately marking the position of the melting curve of helium (T^ — 
38.6 K) at this pressure. 

In Fig. 2, Tc is plotted versus applied pressure to 0.7 GPa and is seen to follow a 
highly linear dependence dT^/dP ~ — 1.11±0.02 K/GPa. The first data point 1' was 
obtained after first applying a pressure of ~ 0.7 GPa at RT before cooling down to low 
temperatures (~ 60 K) and reducing the pressure, yielding Tc ~ 38.88 K at 0.341 GPa. 
Point 2' was measured after releasing the pressure at low temperature, giving Tc ~ 
39.25 K at ambient pressure (0 GPa); no change in occurred after intermittently 
warming the sample to RT (point 3). Further data were obtained following pressure 
changes at both RT (unprimed data) and low temperature (primed data). As is 
observed for the vast majority of superconducting materials without phase change, 
the dependence of Tc on pressure for MgB2 is single-valued and does not depend 
on the pressure/temperature history of the sample; such history effects do occur in 
certain high-T^ oxides containing mobile species at RT||23|. We thus find that for 
He-gas pressure changes at both ambient and low temperature, Tc{P) for MgB2 is a 
linear, reversible function of pressure to 0.7 GPa [p^ . 

In the present experiments the sample is surrounded by fluid helium near Tc ^ 
39 K for all data taken at pressures P < 0.50 GPa so that the slope dTc/dP ~ 
-1.11 K/GPa gives the true hydrostatic pressure dependence of Tc. The fact that 
the sample is in solid helium for P > 0.5 GPa is seen in Fig. 2 to have no effect 
on the pressure dependence Tc{P); indeed, solid helium is the softest solid known. 



Our value of dTc/dP differs significantly from those of other groups[0, |T8|, [T^ (see 
discussion above) using pressure media which are either solid at RT or freeze upon 
cooling down at temperatures well above Tc. 

In view of the strong compressibility anisotropy[|l^] and the sizeable anharmonic- 
ity and non-linear electron-phonon coupling|p] anticipated for MgB2, it is likely that 
shear stresses of sufficient magnitude will cause appreciable changes in the pressure 
dependence of Tc, as observed for other anisotropic substances such as the super- 
conducting oxides p5[ and organic superconductors [^]. The largest shear stresses 
are generated by changing the pressure on a solid pressure medium, such as steatite, 
or using no pressure medium at all. The shear stresses generated in cooling Fluo- 
rinert or other comparable liquids through the melting curve are admittedly much 
smaller and depend on details of the individual experimental procedures used, such 
as the cooling rate, change in applied force upon cooling, etc. Only experiment 
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can determine whether or not the Tc{P) dependences measured in the Fluorinert 
experiments |T7|, [T^ are influenced by shear stresses. To exclude the possibility of 
sample-dependent effects, such experiments should be carried out on a single sample. 
Lorenz et al.^^ have very recently carried out He-gas high-pressure studies on the 



same sample studied by them previously in Fluorinert and find a value of dTc/dP 
equal to their previous value (-1.6 K/GPa), within experimental error; further mea- 
surements in the same He-gas system on a second, high-quality sample yielded the 
dependence dTc/dP ~ —1.07 K/GPa, a value very close to our present result. This 
finding lends support to the observation by Monteverde et al. |^ that dTc/dP in 
MgB2 may be sample dependent. 

The present Tc{P) studies and parallel high- pressure structural studies by Jor- 



gensen et a/.|]T2| were both carried out on the same MgB2 sample over the same 
He-gas pressure range, thus allowing an accurate determination of the change in Tc 
with lattice parameter. The change in Tc with unit cell volume, for example, is given 
by 

using the above values dTjdP ~ -1.11 ± 0.02 K/GPa, B = 147.2 ± 0.7 GPa, and 
Tc = 39.25 K. 

We will now discuss the implications of this result for the nature of the supercon- 



ducting state in MgB2. First consider the McMillan equation p8[ 
Tc ^ {(uj) /1.20) exp {[-1.04(1 + A)] / [A - /i*(l + 0.62A)]} , valid for strong coupling 
(A < 1.5), which connects the value of Tc with the electron-phonon coupling param- 
eter A, an average phonon frequency (uj) , and the Coulomb repulsion //*, which we 
assume to be pressure independent. Taking the logarithmic volume derivative of Tc, 
we obtain the simple relation 

dlnTc . f dlnn ^ ^ 

-7 + A —-L + 27 L (2) 



dlnV ' [dlnV 

where 7 = —din (u) /dlnV is the Griineisen parameter, 77 = N{Ef) (P) is the Hop- 
field parameter [p^] given by the product of the electronic density of states and the av- 
erage squared electronic matrix element, and A = 1.04A [1 + 0.38/i*] / [A — /i*(l + 0.62A)] 

Eq. (2) has a simple interpretation. The first term on the right, which comes 
from the pref actor to the exponent in the McMillan expression for Tc, is usually small 
relative to the second term, as will be demonstrated below. The sign of the loga- 
rithmic derivative dlnTc/ dlnV , therefore, is determined by the relative magnitude 
of the two terms in the curly brackets. The first "electronic" derivative is nega- 
tive {dlnrj/dlnV ~ —1 for simple metals (s,p electrons) [^], but equals -3 to -4 
for transition metals (d electrons) p9|), whereas the second "lattice" term is positive 
(typically 27 3 — 5). Since in simple- metal superconductors, like Al, In, Sn, and 
Pb, the lattice term dominates over the electronic term, and A is always positive. 
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the sign of dlnTc/dlnV is the same as that in the curly brackets, namely positive; 
this accounts for the universal decrease of with pressure due to lattice stiffening in 
simple metals. In selected transition metals the electronic term may become larger 
than the lattice term, in which case Tc would be expected to increase with pressure, 
as observed in experiment ||29||. 

Let us now apply Eq. (2) in more detail to a canonical BCS simple-metal super- 
conductor. In Sn, for example, Tc decreases under pressure at the rate dTc/dP ^ 



-0.482 K/GPa which leads to d\nTjd\nV ~ +7.2 H]. Inserting Tc(0) ~ 3.73 



K, (uj) ~ 110 K pT|, and n* = 0.1 into the above McMillan equation, we obtain 



A ~ 0.69 from which follows that A ~ 2.47. Inserting these values into Eq. (2) and 
setting dlnrj/dlnV ~ —1 from above for simple metals, we can solve Eq. (2) for 
the Griineisen parameter to obtain 7 ~ +2.46, in reasonable agreement with exper- 



iment for Sn (7 ^ +2.1)||20[. Similar results are obtained for other simple-metal 
superconductors. 

We now repeat the same calculation with the McMillan equation for MgB2 using 
the logarithmically averaged phonon energy from inelastic neutron studies]^ (uj) = 
670 K, Tc(0) ~ 39.25 K, and /i* = 0.1, yielding A ~ 0.90 and A ~ 1.75. From Eq. (1) 
we have dlnTc/dlnV = +4.16. Since the pairing electrons in MgB2 are believed to 



be s,p in character ||32||, we set dlwq/dlnV ~ —1 [Q. Inserting these values into Eq. 
(2), we find 7 ~ 2.36, in reasonable agreement with the value 7 ~ 2.9 from Raman 



spectroscopy studies|T^ or 7 2.3 from ab initio electronic structure calculations on 
MgB2 @]. If, on the other hand, one were to use the same bulk modulus but the 
pressure derivative dT^/dP ~ —2.0 K/GPa from Saito et aiWEl, one obtains from 
Eq. (2) the unusually high value 7 ~ 3.7. 

In extensive specific heat0 and high- resolution photoemission studies]^ on MgB2, 
evidence is found for a multicomponent superconducting gap; the latter study also 
reports an inconsistency in the values of the electron-phonon coupling constant from 
McMillan's equation and the renormalization of the electronic density-of-states. These 
results call into question the suitability of the isotropic McMillan equation for describ- 
ing this system. 

Under the above assumptions, we thus conclude that the rate of decrease of Tc with 
pressure found in the present experiments on MgB2 is consistent with BCS phonon- 
mediated superconductivity. The authors hope that the accurate determination of 
the volume dependence of Tc in this work, dlnTc/dlnV = +4.16 + 0.08, will stimulate 
ab initio theoretical calculations. At first glance the present results appear to be 
inconsistent with the model of Hirsch and Marsiglio PH], |TT|] which predicts that Tc 
should increase with pressure for optimally doped samples. However, within their 
model, Tc for a non optimally doped sample may decrease if sufficient change in 
the carrier concentration occurs when pressure is applied. Further studies, such as 
high-pressure Hall effect measurements, are necessary to clarify this possibility. 

Acknowledgments 



6 



Work at Washington University supported by NSF grant DMR 98-03820 and that 
at the Argonne National Laboratory by the U.S. Department of Energy, Office of 
Science, contract No. W-31-109-ENG-38. 



References 

[1] J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, and J. Akimitsu, Nature 
410, 63 (2001). 

[2] J. G. Bednorz and J. A. Miiller, Z. Phys. B 64, 189 (1986). 

[3] S. L. Bud'ko, G. Lapertot, C. Petrovic, C. E. Cunningham, N. Anderson, P. 
C. Canfield, preprint cond-mat/0101463 ; D. G. Hinks, H. Glaus, J. Jorgensen, 



Nature (in press), see preprint in |cond-mat /0104242 



[4] R. K. Kremer, B. J. Gibson, K. Ahn, preprint |cond-mat / 102452 



[5] Y. Wang, T. Plackowski, and A. Junod, preprint |cond-mat/010318T 



[6] T. Takahashi, T. Sato, S. Souma, T. Muranaka, J. Akimitsu, preprint |cond 



mat/0103079 . 



[7] R. Osborn, E. A. Goremychkin, A. I. Kolesnikov, D. G. Hinks, preprint cond- 
mat/010306^ . 



[8] T. Yildirim, O. Giilseren, J. W. Lynn, C. M. Brown, T. J. Udovic, H. Z. Qing, 
N. Rogado, K. A. Regan, M. A. Hayward, J. S. Slusky, T. He, M. K. Haas, P. 
Khalifah, K. Inumaru, R. J. Cava, Phys. Rev. Letters (in press), see preprint 
|cond-mat/0103469| . 



[9] See, for example, P. E. Seiden, Phys. Rev. 179, 458 (1969). 
[10] J. E. Hirsch, preprint |cond- mat / 1 2 1 T5 . 



[11] J. E. Hirsch and F. Marsiglio, preprint |cond-mat/010247^ . 

[12] J. D. Jorgensen, D. G. Hinks, S. Short, Phys. Rev. B (Rapid)(in press). 

[13] I. Loa and K. Syassen, preprint |cond-mat / 102462 . 



[14] A. F. Goncharov, V. V. Struzhkin, E. Gregoryanz, J. Hu, R. J. Hemley, Ho-k. 
Mao, G. Lapertot, S. L. Bud'ko, P. C. Canfield, preprint |cond-mat /0 1 04042 . 



[15] K. Prassides, Y. Iwasa, T. Ito, Dam H. Chi, K. Uehara, E. Nishibori, M. Takata, 



S. Sakata, Y. Ohishi, O. Shimomura, T. Muranaka, J. Akimitsu, preprint cond- 



mat/0102507 . 



7 



T. Vogt, G. Schneider, J. A. Hriljac, G. Yang, J. S. Abell, preprint |cond- 
mat/0102480| . 



B. Lorenz, R. L. Meng, C. W. Chu, Phys. Rev. B (in press); preprint |cond- 
mat/0 102261 . 



E. Saito, T. Taknenobu, T. Ito, Y. Iwasa, K. Prassides, and T. Arima, J. Phys.: 
Condens. Matter 13, L267 (2001). 

M. Monteverde, M. Niinez-Regueiro, N. Rogado, K. A. Regan, M. A. Hayward, 
T. He, S. M. Loureiro, R. J. Cava, Science 292, 75 (2001). 

A. Eihng and J. S. Schilhng, J. Phys. F 11, 623 (1981). 
I. L. Spain and S. Segall, Cryogenics 11, 26 (1971). 

J. S. Schilhng, J. Diederichs, S. Klotz, R. Sieburger, in: Magnetic susceptibility of 
Superconductors and Other Spin Systems, edited by R. A. Hein, T. L. FrancaviUa, 
D. H. Liebenberg (Plenum Press, New York, 1991), p. 107. 

See, for example, S. Sadewasser, J. S. Schilling, and A. M. Hermann, Phys. Rev. 
B 62, 9155 (2000). 

In a separate experiment on MgB2 to 0.077 GPa pressure in Neon gas, an identical 
Tc(P) dependence was obtained. 

S. Klotz and J. S. Schilling, Physica C 209, 499 (1993); J. S. SchiUing and S. 
Klotz, in: Physical Properties of High Temperature Superconductors, Vol. HI, ed. 
D. M. Ginsberg (World Scientific, Singapore, 1992) p. 59. 

T. Ishiguro, H. Ito, Y. Yamauchi, E. Ohmichi, M. Kubota, H. Yamochi, G. Saito, 
M. V. Kartsovnik, M. A. Tanatar, Yu. V. Sushko, G. Yu. Logvenov, Synthetic 
Metals 85, 1471 (1997); D. Jerome and H. J. Schulz, Adv. Phys. 31, 299 (1982). 

B. Lorenz, R. L. Meng, and C. W. Chu, preprint |cond-mat/0104303| . 

W. L. McMillan, Phys. Rev. 167, 331 (1968); P. B. Allen and R. C. Dynes, Phys. 
Rev. B 12, 905 (1975). 

J. J. Hopfield, Physica 55, 41 (1971). 

J. S. Schilling and S. Klotz, in: Physical Properties of High Temperature Super- 
conductors, Vol. HI, ed. D.M. Ginsberg (World Scientific, Singapore, 1992) p. 
59. 

[31] R. C. Dynes, Sohd State Commun. 10, 615 (1972). 



8 



[32] J. Kortus, I. I. Mazin, K. D. Belashchenko, V. P. Antropov, L. L. Boyer, preprint 
|cond-mat/0101446| ; J. M. An and W. E. Pickett, preprint |cond-mat/0102391| ;Y. 



Kong, O. V. Dolgov, O. Jepsen, O. K. Andersen, preprint |cond-mat/010249£ 



N. I. Medvedeva, A. L. Ivanovskii, J. E. Medvedeva, A. J. Freeman, preprint 
|cond-mat/0103157| . 

[33] The volume dependence of the Hopfield parameter dlnrj/dlnV ~ —1 used here 
is in reasonable agreement with that dlnrj/dlnV = Bdlnrj/dP ^ —0.81 ob- 
tained from first-principles electronic structure calculations on MgB2 by N. I. 
Medvedera, A. L. Ivanovskii, J. E. Medvedeva, A. J. Freeman, and D. L. No- 
vokov (preprint |cond-mat/010434^ ) which obtain dlnrj/dP ^ +0.55 %/GPa 



and setting B = 147.2 GPa from Ref. 12. The resulting change in the estimate 
7 ~ 2.23 is small. 

[34] D. Roundy, H. J. Choi, H. Sun, S. G. Louie, M. L. Cohen, post deadline session 
on MgB2, Monday evening, March 12, 2001, APS March Meeting in Seattle, 
Washington. 

[35] S. Tsuda, T. Tokoya, T. Kiss, Y. Takano, K. Tongano, H. Kitou, H. Ihara, and 
S. Shin, preprint |cond-mat / 1 04489 . 



9 



1 Figure Captions 



Fig. 1. Real part of the ac susceptibility of MgB2 versus temperature at ambient and 
high pressures. The applied magnetic field is 0.113 Oe (rms) at 1,023 Hz. Intercept 
of straight tangent lines defines superconducting onset at ambient pressure Tc(0) ~ 
39.25 K. No correction is made for demagnetization effects. 

Fig. 2. Superconducting transition temperature versus applied pressure. Numbers 
give order of measurement. Data for pts. 2', 6, 8', and 11 are shown in Fig. 1. A 
typical error bar for Tc (±0.01 K) is given in lower left corner; the error in pressure is 
less than the symbol size. Pressure was either changed at RT (unprimed numbers) 
or at low temperatures ~ 60 K (primed numbers). 
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